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Abstract 

In this paper we give the first method for constructing n-multimagic squares (and 
hypercubes) for any n. We give an explicit formula in the case of squares and 
an effective existence proof in the higher dimensional case. Finally we prove that 
n-multimagic squares do not exist for certain orders. 

Introduction 

Magic squares have been studied over 4000 years. Recently some exciting new results 
have been found considering these squares. For instance the first method of constructing 
all most-perfect magic squares and their enumeration appeared in (See also [H|). Also 
the non-existence of a 8 x 8 magic knights tour (JSj) has been verified by computers. 
Finally a very natural connection has been found between the properties and construction 
of Franklin's squares and his magic circles ([IDl)- However, there are still many unsolved 
problems (e.g. [l], [12]) and as Clifford Pickover says "the field of magic square study 
is wide open" (|12j. p. 26). In this paper we will concentrate on so called multimagic 
squares. 

Suppose that M is an m x m matrix M consisting of natural numbers. Then M is 
called a magic square if the sum of all elements in each column, row and main diagonal 
gives the same number; the so-called magic number. Let M* d be the matrix obtained by 
raising each element of M to the d-th power. The matrix M is called an n-multimagic 
square (where n is a fixed positive natural number) if M* d is a magic square for d = 
1,2, ... ,n. The matrix M is called normal if its matrix elements consist of the consecutive 
integers 1,2,... , m 2 . Throughout this paper we always consider normal magic squares 
(of course if m > 1 and d > 1 then the matrix M* d is not normal). 

The first 2-multimagic square was published by Pfeffermann in 1891: it has order 
8 (Figure HJ [TT], [5]). In 1905 the first 3-multimagic square was constructed by Tarry: 
it has order 128. In 2001 both a 4- and a 5- multimagic square were constructed by 
Boyer and Viricel respectively of order 512 and 1024 (|3], [3] where they also give a nice 
history on the subject). The record up to now was a 6-multimagic square of order 4096 
constructed by Pan Fengchu in October 2003 ([E])- 
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Figure 1: G. Pfeffermann, 1891 



In this paper we give a constructive procedure to make a large class of n-multimagic 
squares for each positive integer n > 2. The problem of finding such squares is reduced 
to an easy linear algebra problem which is solved in general in section 0] A more explicit 
solution is described in section El This solution is used to give an explicit formula for 
n-multimagic squares for all n > 3. In particular it gives the first 7-multimagic squares, 
of order 13 7 and 8-multimagic squares of order 17 8 etc. 

The method described for constructing n-multimagic squares can easily be extended 
to n-multimagic cubes and hypercubes. We refer to section 0] for all definitions and more 
details. 

1 Preliminaries 

Throughout this paper R denotes a finite ring with q elements. Let R* denote the set of 
elements in R which have a multiplicative inverse. 

Definition 1.1 For c G R we call a bijection N : R, —> {0, 1, . . . , q — 1} of type c if 
N(a) + N(-a + c) = q - 1 for all a G R. 

Lemma 1.2 

i). If 2 G R*, then for every c G R there exists a bijection N of type c. 

ii). If 1 (ji R* , then for every c G R* there exists a bijection N of type c. 

Proof. For c G R define (p = ip c : R — > R by tp(a) = —a + c for all a G R. Then ip 2 is 
equal to the identity. So all orbits of <p have length 1 or 2. An element a is a fixed point 
of <p if and only if 2a = c. For a G R denote its orbit under ip by 0(a). 

i) Let 2 be a unit in R. Then ip has exactly one fixed point, namely oq := 2~ 1 c. The orbit 
O(ao) of do has one element. Let 0(a\), 0(a 2 ), . . . , 0(a s ) be the other orbits. They all 
have two elements. In particular, we get s = (q — l)/2. Define N(a ) = s, iV(oj) = i — 1 
and N((p(a,i)) = q — 1 — N(ai) (= q — i), for i = 1, 2, . . . , s. Now iV is as desired. 

ii) Let 2 not be a unit in R. Then for c G R*, <p = <p c has no fixed points. Indeed, if 
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<p(a) = a, then we get 2a = c G R*. So 2 G if* and we have a contradiction. We can 
write R = IJi=i^( a «)' where each O(aj) is an orbit with two elements. In particular, 
s — q/2. Define iV(dj) = i — 1 and N(ip(a,i)) = q — 1 — iV(ai) (= q — i) for all i. Then iV 
is as desired. □ 

Let m be a positive integer. For each 1 < j < m we choose a bijection 

%):£-{(), l,...,<z-l} 

of type Cj, for some Cj G -R (this is possible by Lemma fOj) . Put c = (ci, . . . , c m ) G -R m 
and define N m : i? m -> {1, 2, . . . , g m } by 

m 

iV TO (ai, . . . ,a m ) = 1 + y^g J ~ 1 iV(j)(aj). 

i=i 

Since the coefficients of the g-adic expansion of any natural number are unique and each 
iVy) is a bijection, it follows that N m is a bijection. 

Lemma 1.3 N m (—a) = q rn + 1 — iV m (a + c) ; for all a = (ai, . . . , a m ) G i? m . 
Proof. From definition 11.11 follows that 

m 

N m {-a) + N m {a + c) = 2 + ^ ^{N^-aj) + N {j) { aj + Cj )) = 

i=i 

= 2 + (1 + g + . . . + g" 1 " 1 )^ - 1) = 2 + (g m - 1) = g m + 1. 

□ 

To conclude this section we will give a result (proposition II. 5|) which plays a crucial role 
in the next section. First some notations. Let n and s be positive integers. Suppose 
that L : R n —>■ R s is an affine map, i.e., there exists an i?-module homomorphism 
Lq : R n — > R s and a vector v G R s such that 

L(a) = L (a) +v, a G -R n . 



Lemma 1.4 7/ L : i? n — > R s is a surjective affine map, then jfL 1 (y) = q n s for all 
y G R s . 

Proof. Let y G R s . Since L is surjective there exists a G R n such that L(a ) = y. 
It follows that L _1 (y) = a + kerL . Since L is surjective, so is L . It follows that 
i? n /kerLo — whence ^kerLo = g ra_s and consequently #L~ l (y) = q n ~ s . □ 
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Proposition 1.5 Suppose that L : R n — ► R s is a surjective affine map. For each j G 
{1,2,..., s}, let N(j\ : R — > {0, 1, . . . , q — 1} be a bijection. For any nonnegative integers 
ei, . . . , e s we have 



%)(^)i) ei • ..N (s) (L(a) s r = q n ~ s ( ) 

aGfi™ \i=0 / \i=0 

Proof. Let y = (yx, . . . , y s ) G R s . Then for each a G we get 

iV^a)!) 61 . . .iV (s) (L(a) s ) e * = iV (1) ( yi ) ei . . . N (s) (y s )^. 
So by Lemma fl.4l we obtain 

N {1) (L(a) i r...N (s) (L(a) s r = q^N^y^ . . . N (s) (y s r . (1) 

Since L is surjective, _R n is the disjoint union of the fibres L _1 (?/), y G R s . We deduce 
from Q that 

'£N {l) (L(a) 1 r...N {s) (L(a) s y° 

a£R n 

= E E N (1) (L(a) i r...N (s) (L(a) s r 



E^ S %)(2/i) ei ---^)(^)' 

y£R° 

We^^Y-Ye 

\yieR / \y s eR 



N is) (y s ) ea 



\yieR / \y a £R 

Since each Nq) : R — > {0, 1, . . . , q — 1} is a bijection we get that 

E%)(^ = X>, 

i=0 

which concludes the proof. □ 



2 A construction of n-multimagic squares 

Let n G N. The following theorem gives the main tool for constructing n-multimagic 
squares. We will use the notations introduced in the previous section. Let R be a finite 
ring with q elements. Write Gl m (R) for the ring of m x m invertable matrices over R. 
First we choose c±, . . . , c n in R and n bijections 

N (1) ,...,N (n) :R^{0,1,. ..,<?-!} 
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of types ci, . . . , c n respectively. With these (not necessarily different) bijections we define 

N n :R n ^{l,2,...,q n } 

as described in section ^ This choice will be fixed throughout this section. In a similar 
way we fix a bijection 

N^ n :R 2n ^{l,2,...,q 2n } 
using bijections -Wm, • • • , NL n -\ of type c[, . . . , c' 2n respectively. 

Definition 2.1 A matrix X G Gl2 n (R) is called an n-multimagic generator matrix if for 
the matrices A,Be Mat2 nn (i?) suc/i i/ia£ X = (A B), all n x n minors of A, B, A + B 
and A — B are units in R. 

Theorem 2.2 Suppose that X G Gl2 n {R) is an n-multimagic generator matrix. For any 
t G R 2n the q n x q n matrix M defined by 

M Nn(a) , Nn(b) =NL(xl a ')+t), a,beR n 



is n-multimagic. 



Note that the matrix M in the theorem is well-defined because N n (a) takes each value 
in {1,2, ... ,q 2 } exactly once, and so does N n (b). 

Proof. 

i) First observe that all matrix elements My are distinct since X e Gl2 n {R) and N' 2n 
is a bijection. Consequently the matrix M consists exactly of all elements of the set 
{1, 2, . . . , q 2n }, so M is normal. 

ii) Now let d be an integer with 1 < d < n and write X = (A B) with A,B& Mat2 n ,n(-R)- 
First we want to show that the sum of all elements in any column of M* d is the same 
constant which only depends on q and n. Fix b G R n . Then the N n (b)-t\\ column of 
M* d consists of the elements ^ N ^, where a runs through R n (remember that 
N n : R n -> {1, 2, . . . , q n } is a bijection)." 

Let Sb(d) be the sum of the elements of the iV n (&)-th column of M* d , so 

S b (d) = M N n (a),N n (b)- 
a£R n 

To compute Sb(d) first observe that the j-th component of the vector X{Vj = (A B)(V) 
is equal to A^ ■ a + B^ ■ b, where A(j) (respectively B^) denotes the j-th row of A 
(respectively B). Using the definitions of M Nn ^ a )^ n {b) and we get 

(2n \ d 

l + J^Cfab)] (2) 
i=i / 
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where 

Cj(a, b) = g- ? " 1 iV ( ' j) (A (j) .a + B {j) .b + tj), for all 1 < j < 2n. (3) 

Now observe that + - • ■ + X2 n ) d can be written as 1+g, where g is a sum of terms of 

the form ax^ . . . Xj* , where 1 < j\ < ji ■ ■ ■ < js < 2n, t\ , . . . , e s > 1 and e\ + • ■ ■ + e s < d 
(so in particular s < d < n) and a is a positive integer. So it follows from (J2J that Sb{d) 
only depends on q and n if we can show that for each set of exponents e±, . . . ,e s and 
indices ji, ■ ■ ■ ,j s as above, the sum 

J2C h (a,br...C js (a,by° (4) 

aeR n 

only depends on q and n (and of course e%, . . . , e s , ji, . . . , j s ). To see this we are going to 
use Proposition ll.5l Therefore put J = (Ji, ■ ■ ■ ,j s ) an d define the affine map L : R n — > R s 
by the formula 

L(a) = • a + B {J) -b + tj 

where (respectively -B(j)) is the s x n matrix with rows Ay^, . . . , A(j s ) (respectively 
Bfa), . . . ,B(j a )) and in the column of length s with components tj 17 . . .tj g . Since, 
as observed above, s < n and all n x n minors of A are units in R (by hypothesis) it 
follows that L : R n — > R s is surjective. By (0) we have C^(a, 6) = q ji ~ 1 N'^^(L(a) i ) for 
all 1 < z < s. It follows from Proposition II .51 that the expression in (J3J) is equal to 

q n ~ S ■ g ei ^ 1_1 ) H l-esO's- 1 ) . 

which indeed only depends on g and n, as desired. 

iii) Interchanging the roles of a and b in the argument given in ii) we get that all rowsums 
of M* d are equal to the same constant. 

iv) Now let us compute the sum of the (main) diagonal elements of M* d . This sum is 
equal to 

D(d) = J2 K n (a),N n (ay 
a&R n 

To compute D(d) we just repeat the arguments given in ii) with b = a. It then remains 
to show that the expression in (jlj) with b = a equals the expression given in (jSJ), since 
this results in D(d) being equal to the same constant as the rows. Therefore just observe 
that Cj^a, a) = q^' l N' ( jALi{a)i) for all 1 < i < s, where L\ : R n — > R s is the affine 
map defined by 

L(a) = A(j) ■ a + B(j) ■ a + t (J) = (A + S) (J) • a + t (J) 

(recall that J = (ji, ■ ■ ■ ,j s ) )• Since by hypothesis all n x n minors of A + B are units in 
R, it follows that L is surjective. Then using Proposition 11.51 again we obtain that the 
expression in ^ with b = a is indeed equal to the expression given in ©. 



&)■■■(¥) 



6 



v) Finally we compute the sum of all elements from the "second" diagonal of M* d . This 
sum D\d) is equal to 

D '(d) = 2J M N n (a),q™+l-N n (a)- 
a£R n 

Since by lemma fOl q n + 1 — N n (a) = N n (—a + c) we get 

D\d) = J2 M d Nn{alNn( . a+c) . 

a£R™ 

Then repeating the arguments in ii) with b replaced by —a + c leads us to define the 
affine map L' : R n — > R s by 

L'(a) = A (J) ■ a + B {J) (-a + c) + t {J) = (A - B) {J) ■ a + (B (J) ■ c + t (J) ). 

The map V is surjective since all n x n minors of A — B are units in R. So again we 
find that the expression in Q with b replaced by —a + c is equal to the expression in 
(0), resulting in D'(d) being equal to the same constant, which completes the proof of 
this theorem. □ 



3 Finding Generator Matrices 

In order to be able to construct effectively n-multimagic squares by the method described 
in theorem !2.2| we need to show how to find a ring R and a n-multimagic generator matrix 
X G Gl2n{R) which satisfy the conditions of that theorem. 

Below we describe an explicit construction of (n-multimagic) generator matrices. For 
an effective existence proof in a more general setting see section 14.21 

Lemma 3.1 Let n > 1 and R a ring such that 2 and 3 are units in R. If A is an2nxn 
matrix such that every nxn minor of A is a unit in R, then there exists an2nxn matrix 
B such that (A B) is an n-multimagic generator matrix. More precisely, if A = (Q with 
P,Q G M&t n ^ n (R), then we can take B = (_2q)- 

Proof. Since 2 is a unit in R the hypothesis on P and Q implies that also the nxn 
minors of B are units in R. Furthermore the nxn minors of A + B = (^) are also 
units in R, since 3 and —1 are. Similarly the nxn minors of A — B are units in R. 
Finally, using elementary column operators one can reduce the matrix (A B) = (g _^g) 

to the matrix _ 4 °) which is clearly invertible over R since both detP and det(— 4Q) 
are units in R. □ 

Lemma 3.2 Let n > 2 and define the 2n x n matrix A by Aij = (i — for i = 

1, 2, . . . , 2n - 1 and j = 1,2,..., n (0° = 1), A 2n ,j = for j'= 1, 2, . . . , 2n - 1 and 
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.4 



2n,n 



1. So we have 



A 



( 1 
1 
1 
1 




1 

2 
3 




1 

4 
9 



V 



(2n - 2) (2n - 2) 2 





1 

2 n-l 
3 n-l 

(2n - 2)™- 1 

1 



/ 



7/i2 a nn^ i/ien we may mew A as a matrix with entries in R. If{l, 2, . . . , 2n— 2} C 
t/ien all n x n minors of A are units in R. 

Proof. Using Vandermonde determinants one easily verifies that each factor appearing 
in each n x n minor of A is of the form i — j where 0<j<i<2n — 2, from which the 
desired result follows. □ 

As an immediate consequence of the lemmas above we get 

Corollary 3.3 Let n > 2 and R be a ring such that {3, 1, 2, . . . , 2n — 2} C R* , then the 
matrix (A B) is an n-multimagic generator matrix where A is as in Lemma \S. S\ and B 
is as in Lemma \H. 11 



In theorem 13.41 below we choose one bijection iV : R — > {0, 1, 
c G R and define for each m > 1 



q — 1} of some type 



iV m (oi, . . . , dm) = 1 + ^ lN ( a j)- 



3=1 

So in the definition of N m as given in section ^ we take all to be equal to N. 
Theorem 3.4 (An explicit formula for n-multimagic squares.) 

Let n > 3, q a prime number > 2n — 1, R — W q and N : R — > {0, 1, . . . , q — 1} the 
bijection (of type —1) given by N(i) = i for all < i < q — 1. Le£ A and i? 6e as in 
\H. 'A Then for X = (A B) and each t G R 2n the matrix M defined in \2.%A is n-multimagic. 



In other words, for every positive integer n there exists a (normal) n-multimagic square. 
In particular for n = 7 we get 7-multimagic squares of orders 13 7 , 17 7 , 19 7 etcetera. For 
n = 8we get 8-multimagic squares of orders 17 s , 19 s , 23 8 etcetera. 



4 Multimagic cubes and hypercubes 

In this section we briefly indicate how the method developed in the previous sections can 
be extended to construct multimagic cubes, perfect multimagic cubes and hypercubes. 
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4.1 Multimagic cubes 



Note that there is no consensus on the definition of multimagic cubes (hypercubes etc.) 
in the literature. The choice given below can also be found in ^H], 

A cube of numbers (respectively the consecutive numbers 1,2,..., n 3 ) is called magic 
(respectively normal magic) if the sum of all elements in each row, column and pillar 
is the same and is equal to the sum of all elements of each of the four space diagonals. 
Furthermore, if n > 1 such a cube is called n-multimagic if for each 1 < d < n the cube 
obtained by raising each of its elements to the <i-th power is magic. 

Completely analogues to the construction of n-multimagic squares in 12.21 we define a 
q n x q n x q n cube by the formula 



where each of the vectors a, b and c runs through R n , t e R 3n and A, B and C are matrices 
in Mat3„ in (i2) which satisfy the following properties (which guarantee the matrix M to 
be an n-multimagic cube): 

1. (A B C) G Glsn{R) (which guarantees that all the natural number 1, 2, ... , q 3n (= 
(q n ) 3 ) appear in M). 

2. all n x n minors of the matrices A, B and C are units in R (which guarantees that 
for each 1 < d < n the sum of all elements in each column, row and pillar of M* d 
is the same, and hence equal to the magic sum). 

3. all n x n minors of the matrices A + B + C, -A + B + C, A-B + C and A + B-C 
are units in R (which guarantees that for each 1 < d < n the sum of all elements 
on each of the four space diagonals of M* d is equal to the magic sum). 

Recall that a magic cube is called perfect if additionally the diagonals of each orthogonal 
slice have the magic sum property. Furthermore, if n > 1 such a cube is called n- 
multimagic perfect if for each 1 < d < n the cube obtained by raising each of its 
elements to the d-th power is perfect magic. 

To guarantee that a n-multimagic cube M as above is also n-multimagic perfect we 
impose on the matrices A, B, C the following conditions 

4. all n x n minors of the matrices A + B, A — B, A + C, A — C, B + C and B — C 
are units in R. 

4.2 More Generator Matrices 

To find a ring R and matrices A, B and C satisfying the properties 1, 2, 3 and 4 one can 
use the method described below. To facilitate generalizations to (even) higher dimensions 
we give a more general notion of generator matrix: 



M Nn{a) , Nn{b) , Nnic) = N' 3n \(A B CO b \+t 
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Definition 4.1 Let d be a positive integer > 2. We call a matrix X e Gld n (R) an 
n-multimagic d-generator matrix if when we write X = (A 1 . . . A d ) with Ai, . . . , A d e 
Matrf ni „(-R) 7 we have that for 81, . . . , e { — 1, 0, 1} not all 0, the matrix 



has all n x n minors in R* . 

Note that the four properties of the matrix (A B C) of the previous section are equivalent 
to saying that it is a n-multimagic 3-generator matrix. 

Example 4.2 A 1-multimagic 3-generator matrix with R = ¥ q with q prime and q > 11 




Example 4.3 A 2-multimagic 3-generator matrix with R = ¥ q with q prime and q > 11 



(\ 
1 1 
1 
1 
1 



2 
2 
2 
2 
-2 





2 
4 
6 



4 
4 

-4 

-4 



\ 

4 
-8 
-12 



-8 -4 -16 

-2 o -4 y 



Example 4.4 A 2-multimagic 2-generator matrix with R = ¥ q with q prime and q > 
11, q ^ 17: 

/ 1 -1 -2 -1 \ 
10-1-1-2 
1-2-1-1 

12 1 1 

112 1 

\ 2 1 1 1 / 

The following (well-known) lemma is the crucial tool in finding effectively ci-generator 
matrices for all d. 



Lemma 4.5 Let m > 1 and Q(xi,...,x m ) a non-zero polynomial in the variables 
Xi,...,x m overZ. Then one can effectively find integers a±, . . . ,a m such thatQ(ai, . . . ,a m ) 

7^0. 

Proof. We will use induction on m. The case m — 1 is obvious since Q{x\) has at most 
deg Q zeros. Now let m > 2. Write Q = qd(xi, • • • , x n -\)x d m + ■ — \-qi(xi, . . . , x m -\)x m + 
qo(x\, . . . , x m -i) with q d ^ 0. By the induction hypothesis there exist integers a±, . . . , a m -i 
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such that g<i(ai, . . . , a m _i) 7^ 0. So the polynomial q(x m ) = Q(ai, . . . , a m _i,x m ) in Z[x m ] 
is non-zero and has x m -degree d. Consequently there exists an integer a m such that 

q(a m ) = Q(ai, . • • , a m _i, a m ) 7^ 0, 

as desired. □ 

Proposition 4.6 Lei n > 1 and d > 2. Then one can effectively find a positive in- 
teger q > 1 and matrices A\, . . . ,Aa G Matd„ in (Z/gZ) swc/i t/ia£ X = (Ai . . . A4) G 
G/rf n (Z/gZ) an n-multimagic d-generator matrix. 

Proof. To avoid complicating an easy matter we only give a proof for d = 2. For other 
d the procedure is similar. 

Let A u = (Aij) and B u = (Bij) be two universal 2nxn matrices, i.e., the entries Aij and 
Bij are distinct variables. Then each nxn minor of A u , B u , A u +B u and A u — B u is a non- 
zero polynomial in the 4n 2 - variable polynomial ring %[Aij,Bij, 1 < % < 2n, 1 < j < n]. 
Let P be the product of all these minors and let Q be the product of P and the polynomial 
det(A u B u ). By lemma H~B1 we can find integers cijj and b^j such that Q(aij,b it j) is a 
non-zero integer. Finally let q be a positive integer > 1 such that gcd(Q(a it j, g) = 1 
for all i,j. Then one easily verifies that X = (a^) and B = (bij) represent matrices in 
Mat2n,n(Z/gZ) having the desired properties. □ 

4.3 Perfect Multimagic Hypercubes 

From the above it is now clear how to generalize these definitions and constructions to 
higher dimensional hypercubes. 

Theorem 4.7 For all integers d > 2 and n > 1 there exists a normal perfect n- 
multimagic d- dimensional hypercube. 



5 Orders 

We take a short look at the possible orders a n-multimagic square might have. First we 
will show that we can use n-multimagic squares to construct new n-multimagic squares 
of different orders. 

Definition 5.1 Let A G Mat mjTO (i2) and B G Mat riin (i2). We define A-kB G Mat mn>TOn (i2) 
by 

A * B m k+i t mi+j = rn B^i + A^ 
where < i, j < m — 1 and < k, I < n — 1 . 

Proposition 5.2 If A G Mat m>m (i2) and B G Mat ri) „(i2) are p -multimagic squares then 
A-k B is an p-multimagic square. 
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Proof. The proof is a straightforward calculation so we only explicitly show that the 
columnsum is independent of the columnindex. We use the same notations as in the 
definition above. Fix 1 < e < p. We denote S x = J2i^i,j an d T x = Yli^fj f° r 
1 < x < p. Since A and B are p-multimagic squares these S x and T x are constants. Fix 
the columnindex (5 — ml + j and write a = mk + i. 

mn—l m—1 rt— 1 

E(^ B )^ = EE( ffl ^ + ^) e 

a=0 fe=0 i=0 

= EEE(f) ^ = E E (°) ™ 2, "^r E * 
= E E C) m"<~)«sr = E (I) '" 2, " ,s - E 
= E(l)™ 2<e "" s - r - 

So the columnsum is independent of hence a constant. □ 

From the foregoing one might get the impression that n-multimagic squares are con- 
structable for all (large enough) orders. This is not the case as we will prove next. 

Definition 5.3 If p is a prime and n is a positive integer, then we define v p (n) as the 
largest nonnegative integer e such that p e divides n. For a non-zero rational number ^ 
we define v p (^) = v p (n) - v p (m). 

Lemma 5.4 (basic properties of v p ) 
Let a, e, m, p G N with p prime. 

i) v p (r ■ s) = v p (r) + v p (s) for any non-zero rational numbers r and s 

ii) If v p (m) = e and 1 < a < p e then v p (m 2 — a) = v p (a). 

Proof. Straightforward using the unique factorization of integers. □ 

Theorem 5.5 Let M be a normal n-multimagic square of order m and p GN a prime. 
Ifv p (m) — e > 1 then n < p e+1 — 2. 

Proof. Write M = (m tJ ) !J=1 .. m . Note that if / e C[x], /(Z) C Z and deg x (f) < n 
then the matrix M°f = (/(^jj))jj=i.. m € Mat mjm (Z) is also a (general) magic square. 
Take f(x) = ( x ~ 1 ) £ Q[x] and note that it is a polynomial of degree n with /(Z) C Z. 
Since M°f is magic the sum of all m 2 elements of M°f is an integer which is m times 



12 



the magic sum, i.e. the sum of any row or column. Since /(Z) C Z the magic sum is an 
integer which implies that m divides 



2 2 

x=l x=l ^ 



x — 1\ / m? 
n J V n + 1 



This means that the (rational) number 

m(m 2 — l)(m 2 — 2) • • • (m 2 — n) 
(n + 1)! 



(6) 



is really an integer. 

Since an n-multimagic square of order m can not exist if an (n — l)-multimagic square 
of order m does not exist, it is enough to show that the case n = p e+l — 1 is impossible 
by showing that (P) is not an integer. 

So assume n = p e+1 — 1. Note that © is not equal to since p e divides m, so in particular 
m 2 > p 2e > p e+1 — 1 = n (remember e is positive). Using lemma 15^1 we see that 



v p 



m 



m(m 2 — l)(m 2 — 2) • • • (m 2 — n) 
(n + 1)! 

'm 2 — 1) (m 2 — 2) (m 2 — n) 



n+1 1 2 n 



(use lemma IB~H i) 



( m \ ^> / m 2 — i\ 
+ / v v\ = ( use lemma l5T4l ii) 

rt 



i=i 

rn 



P \n+1 



v 



i=i 
m 



e-(e + l) = -1 

So (JHJ) is not a integer if n = p e+1 — 1. □ 
In particular this shows that there are no 3-multimagic squares of order m = 2 (mod 4). 



6 More Examples 

To conclude this paper we give some interesting multimagic squares using theorem 12.21 
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In 16.11 - 16.31 below we choose one bijection N 
c G R and define for each m > 1 



R 



{0, 1, . . . , q — 1} of some type 



iV m (oi, . . . , a m ) = 1 + Q j ^ N ( a j)- 

i=i 

So in the definition of N m as given in section we take all to be equal to N. The 
notations are as in theorem 12.21 where we take N' 2n = N 2n . 

Example 6.1 (A family of associative bimagic squares of order 16.) 
Take R = ¥ 2 [x]/(x 2 +x + l), 



X 



( x 1 x \ 
111 
1110 
V x 1 x J 



t G -R 4 arbitrary and N : R — > {0, 1, 2, 3} (a bijection of type x + 1) given by iV(0) = 0, 
JV(1) = 2, N(x) = 1 and JV(x + 1) = 3. 

Then the corresponding matrix M defined in theorem 12.21 is bimagic (= 2-multimagic) 
and associative (= the sum of any pair of matrix elements which are symmetric with 
respect to the center of the square is equal to 16 2 + 1). 



41 


252 


74 


155 


125 


176 


30 


207 


129 


84 


226 


51 


213 


8 


182 


103 


62 


239 


93 


144 


106 


187 


9 


220 


150 


71 


245 


40 


194 


19 


161 


116 


3 


210 


100 


177 


87 


134 


56 


229 


171 


122 


204 


25 


255 


16 


160 


77 


24 


197 


119 


166 


68 


145 


35 


242 


192 


109 


223 


14 


236 


57 


139 


90 


240 


61 


143 


94 


188 


105 


219 


10 


72 


149 


39 


246 


20 


193 


115 


162 


251 


12 


156 


73 


175 


126 


208 


29 


83 


130 


52 


225 


7 


214 


104 


181 


198 


23 


165 


120 


146 


67 


211 


36 


110 


191 


13 


224 


58 


235 


89 


140 


209 


4 


178 


99 


133 


88 


230 


55 


121 


172 


26 


203 


45 


256 


78 


159 


98 


179 


1 


212 


54 


231 


85 


136 


202 


27 


169 


124 


158 


79 


253 


48 


117 


168 


22 


199 


33 


244 


66 


147 


221 


16 


190 


111 


137 


92 


234 


59 


76 


153 


43 


250 


32 


205 


127 


174 


228 


49 


131 


82 


184 


101 


215 


6 


95 


142 


64 


237 


11 


218 


108 


185 


247 


38 


152 


69 


163 


114 


196 


17 


167 


118 


200 


21 


243 


34 


148 


65 


15 


222 


112 


189 


91 


138 


60 


233 


180 


97 


211 


2 


232 


53 


135 


86 


28 


201 


123 


170 


80 


157 


47 


254 


141 


96 


238 


63 


217 


12 


186 


107 


37 


248 


70 


151 


113 


164 


18 


195 


154 


75 


249 


44 


206 


31 


173 


128 


50 


227 


81 


132 


102 


183 


5 


216 



Example 6.2 (A family of bimagic squares of odd order. ) 
Take R = Z/gZ with q > 3, q odd, 



X 



( 1 1 \ 

2 1 

112 1 

\ 2 1 2 2 / 



t E R A arbitrary and N : R — > {0, 1, . . . , q — 1} the standard bijection (of type —1) as in 
13.41 Then the corresponding q 2 x q 2 matrices M as defined in theorem 12.21 are bimagic. 
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In particular this gives a family of bimagic squares of odd order. 



1 35 60 
70 14 39 
49 74 27 


23 48 79 
56 9 31 
44 69 10 


18 40 65 
78 19 53 
30 61 5 


38 72 13 
26 51 73 
59 3 34 


33 55 8 
12 43 68 
81 22 47 


52 77 21 
4 29 63 
64 17 42 


75 25 50 
36 58 2 
15 37 71 


67 11 45 
46 80 24 
7 32 57 


62 6 28 
41 66 16 
20 54 76 



Finally, if we choose q = 3 and t=(2120)*we recover the associative 9x9 bimagic 
square constructed by R.V. Heath (see p. 212 PJ) from before 1974. 

Example 6.3 (An associative, pandiagonal, bimagic, .... magic square of order 25.) 
Take R = F 5 , N : R — > {0, 1, 2, 3, 4} the standard bijection (of type —1) and 



X 





1 


2 


2 


1 


2 


2 


4 


1 


3 


3 


4 


1 


4 


3 


2 



4 


V 2 / 



Then the corresponding 25 x 25 matrix is associative, pandiagonal, bimagic and has the 
following properties: 

i) Each of the 25 standard 5x5 submatrices is pandiagonal (even with the same 
magic sum). 

ii) For each pair (i, j) (1 <i,j< 25) the 5x5 matrix obtained by deleting each row 
with row number not equivalent to i mod 5 and each column with column number 
not equivalent to j mod 5 is pandiagonal! 



103 350 567 164 381 
167 389 106 328 575 
331 553 175 392 114 
400 117 339 556 153 
564 156 378 125 342 


291 513 235 452 74 
460 52 299 516 238 
524 241 463 60 277 
63 285 502 249 466 
227 474 66 288 510 


584 176 423 20 362 
23 370 587 184 401 
187 409 1 373 595 
351 598 195 412 9 
420 12 359 576 198 


147 494 86 308 530 
311 533 130 497 94 
480 97 319 536 133 
544 136 483 80 322 
83 305 547 144 486 


440 32 254 621 218 
604 221 443 40 257 
43 265 607 204 446 
207 429 46 268 615 
271 618 215 432 29 


134 476 98 320 537 
323 545 137 484 76 
487 84 301 548 145 
526 148 495 87 309 
95 312 534 126 498 


447 44 261 608 205 
611 208 430 47 269 
30 272 619 211 433 
219 436 33 255 622 
258 605 222 444 36 


115 332 554 171 393 
154 396 118 340 557 
343 565 157 379 121 
382 104 346 568 165 
571 168 390 107 329 


278 525 242 464 56 
467 64 281 503 250 
506 228 475 67 289 
75 292 514 231 453 
239 456 53 300 517 


591 188 410 2 374 
10 352 599 191 413 
199 416 13 360 577 
363 585 177 424 16 
402 24 366 588 185 


290 507 229 471 68 
454 71 293 515 232 
518 240 457 54 296 
57 279 521 243 465 
246 468 65 282 504 


578 200 417 14 356 
17 364 581 178 425 
181 403 25 367 589 
375 592 189 406 3 
414 6 353 600 192 


141 488 85 302 549 
310 527 149 491 88 
499 91 313 535 127 
538 135 477 99 316 
77 324 541 138 485 


434 26 273 620 212 
623 220 437 34 251 
37 259 601 223 445 
201 448 45 262 609 
270 612 209 426 48 


122 344 561 158 380 
161 383 105 347 569 
330 572 169 386 108 
394 111 333 555 172 
558 155 397 119 336 


441 38 260 602 224 
610 202 449 41 263 
49 266 613 210 427 
213 435 27 274 616 
252 624 216 438 35 


109 326 573 170 387 
173 395 112 334 551 
337 559 151 398 120 
376 123 345 562 159 
570 162 384 101 348 


297 519 236 458 55 
461 58 280 522 244 
505 247 469 61 283 
69 286 508 230 472 
233 455 72 294 511 


590 182 404 21 368 
4 371 593 190 407 
193 415 7 354 596 
357 579 196 418 15 
421 18 365 582 179 


128 500 92 314 531 
317 539 131 478 100 
481 78 325 542 139 
550 142 489 81 303 
89 306 528 150 492 


597 194 411 8 355 
11 358 580 197 419 
180 422 19 361 583 
369 586 183 405 22 
408 5 372 594 186 


140 482 79 321 543 
304 546 143 490 82 
493 90 307 529 146 
532 129 496 93 315 
96 318 540 132 479 


428 50 267 614 206 
617 214 431 28 275 
31 253 625 217 439 
225 442 39 256 603 
264 606 203 450 42 


116 338 560 152 399 
160 377 124 341 563 
349 566 163 385 102 
388 110 327 574 166 
552 174 391 113 335 


284 501 248 470 62 
473 70 287 509 226 
512 234 451 73 295 
51 298 520 237 459 
245 462 59 276 523 
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More research into different properties and various examples can be found in the thesis 
of the second author ([HI)- The reader is also referred to the website ([Zj). 



Acknowledgment: The authors like to thank Michiel de Bondt for various stimulating 
and helpful discussions. 
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